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While the interaction potential between two dipoles residing in a single plane is repulsive, in 
a system of two vertically adjacent layers of dipoles it changes from repulsive interaction in the 
long range to attractive interaction in the short range. Here we show that for dipolar excitons in 
semiconductor heterostructures, such a potential may give rise to bound states if two such excitons 
are excited in two separate layers, leading to the formation of vertically coupled dipolar exciton 
molecules. Our calculations prove the existence of such bound states and predict their binding 
energy as a function of the layers separation as well as their thermal distributions. We show that 
these molecules should be observed in realistic systems such as semiconductor coupled quantum well 
structures and the more recent van-der-Waals bound heterostructures. Formation of such molecules 
can lead to new effects such as a collective dipolar drag between layers and new forms of multi¬ 
particle correlations, as well as to the study of dipolar molecular dynamics in a controlled system. 


Two-dimensional (2D) fluids of quantum dipoles 
are currently central to recent theoretical and ex¬ 
perimental efforts to explore many-body physics 
in the quantum regime. Research on such fluids 
ranges from atomic and molecular systems [1, 2] 
to condensed-matter semiconducting heterostruc¬ 
ture systems [3-7]. In the latter, dipolar excitons, 
also known as indirect excitons (IXs), are quasi¬ 
particles composed of a Coulomb-bound electron- 
hole pair, where each of the constituents resides in 
a separate material layer. Excitation of many such 
IXs results in the formation of a two-dimensional 
(2D) layer of oriented dipoles, with a repulsive in¬ 
teraction between particles. Due to their bosonic 
character and relatively long range interaction, these 
particles display a rich phase diagram, with collec¬ 
tive phases and correlation regimes, such as liquid¬ 
ity, dark condensation, and complex spin textures 
[8-14] which are absent in weakly interacting Bose 
gases. In the last decade, there has been a signifi¬ 
cant progress in the formation, control, and under¬ 
standing of dipolar IX quantum fluids in GaAs-based 
bilayers, consisting of biased double quantum well 
(DQW) structures. Recently, new material systems 
such as GaN-based polar quantum wells [15] and 
van-der-Waals (vdW) heterostuctures [16, 17] enable 
formation of dipolar IXs at elevated temperatures, 
with dipole sizes ranging from few to tens of nano¬ 
meters. With such systems, the collective quantum 
phenomena might even be observed up to room tem¬ 
perature [18]. In all these systems the dominant in¬ 
plane repulsion limits the formation of more complex 
structures, such as bi-excitons or charged complexes, 
to a narrow, hard-to-access parameters range and to 
very low binding energies [19-22], suppressing po¬ 


tential observation of interesting phases that may 
arise from such pair correlations. 
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FIG. 1. (a) A system of two dipoles of size di,d 2 resid¬ 

ing in different layers. L is the separation between the 
centers of the dipoles and r is the in-plane distance, (b) 
An illustration of a vertical IXs molecule (dashed oval) 
together with unbound IXs in the two separated layers. 

Here, we propose a double bilayer structure, com¬ 
posed of two vertically stacked DQWs separated by 
a small vertical distance L (see Fig. 1(a)). This 
proximity results in an attractive part for the out 
of plane dipole-dipole interaction between the dipo¬ 
lar IXs formed in each of the bilayers. As we show 
hereafter, this attractive interaction results in the 












formation of vertical IXs molecules, which are com¬ 
posed of bound pairs of two IXs from the two verti¬ 
cally coupled DQWs as illustrated in Fig. 1(b). This 
finding opens new possibilities, such as observing a 
transition from unbound, single IXs gas to a molec¬ 
ular gas when the temperature goes down or as the 
density is modified, and the possibility of attractive 
dipolar drag between IXs in the two coupled layers. 

In the case of two oriented IXs residing in a sin¬ 
gle layer, there is a repulsive force between them 
at all lateral separation distances r. This repulsion 
gives rise to the well-known blue-shift of the exciton 
emission spectrum [9]. However, the electrostatic 
interaction between two oriented IXs residing in dif¬ 
ferent yet proximate layers is qualitatively different: 
for small lateral distance r the two IXs experience 
electrostatic attraction, that gradually changes to 
repulsion for large lateral distances. This Coulomb 
interaction potential Ud{r) between the two IXs is 
given by 
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where di ^2 are the vertical sizes of each of the two 
IXs, L is the vertical separation between the centers 
of the dipoles and 5 is the dielectric constant (it is 
implied here that L > {di + d2)/2). 

The potential in Eq. (1) is plotted in Fig. 2(a) and 
has the following generic features: it reaches a mini¬ 
mum at zero lateral separation where Ud{0) < 0, and 
remains negative up to a distance of the order of L, 
signifying attraction for r ^ L. As for larger dis¬ 
tances where r ^ L, the potential Ud{r) turns repul¬ 
sive and decays as Ud{r) ~ d\d 2 /r‘^ as expected for 
the interaction of the two in-plane aligned dipoles. 
Finally, we note that the potential in Eq. (1) is neu¬ 
tral on average, f d‘^rUd{r) = 0. 

Any net attraction in two dimensions results in 
a bound state [23]. However, the neutral poten¬ 
tial Ud{r) is also known to produce a weakly bound 
state [24, 25]. The question is therefore can experi¬ 
mentally significant bound states of this type exist? 
We thereafter consider two realistic systems of ver¬ 
tically coupled indirect excitons: electrically biased 



r/L 


(b) 



y/L x/L 


FIG. 2. (a) Interaction energy of the two dipoles shown 

in Fig. 1, in units of e^/eL, as a function of the lat¬ 
eral distance measured in units of the vertical separa¬ 
tion between the two quantum double wells, r/L. Here 
the sizes of the dipoles are di/L = 4/9, d 2 /L = 7/18. 
The inset shows a zoom around zero energy to empha¬ 
size the transition from attraction to repulsion beyond 
r ~ L, together with a comparison to a pure dipole- 
dipole repulsive interaction ~ dxd^/r^ (red, dash-dot), 
(b) An exemplary numerical solution of the Schrodinger 
equation showing the probability density of the relative 
in-plane displacement p = {x/L, y/L) of the two exci¬ 
tons in the bound state, with the potential given in (a), 
assuming a GaAs material system with L = 300A. 


GaAs coupled quantum wells and vdW heterostruc¬ 
tures. As we demonstrate below, the binding ener¬ 
gies of these vertically coupled IXs molecules, which 
we denote by IXm, are significant and experimen¬ 
tally relevant. We further argue that the formation 
of complexes which consist of more than two IXs is 
impossible in the considered double bilayer system. 

The Schrodinger equation for the relative in-plane 
center-of-mass separation between the two IXs, p = 
r/L, reads 
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where /i is the reduced mass of the two IXs and Ud 
is the potential given in Eq. (1) with jsL factored 
out. First we estimate the typical scales of the prob¬ 
lem. For GaAs heterostructures L ^ 30nm, while for 
vdW material systems the typical inter-layer sep¬ 
aration is L ^ Inm [26]. The Coulomb interac¬ 
tion is characterized by the energy scale e^jeL ^ 
12meV/(I//10nm) for GaAs-based bilayers (using 
5 « 12), and by e^je'L ^ 230meV/(I//lnm) in 
vdW coupled heterostructures (taking the effective 
dielectric constant s' ~ {1 esi)/2 ^ 6.35 for a 
system placed on a Si substrate with esi = 11.7 
[17]). The kinetic energy scale in GaAs-based bi¬ 
layers is Y? j ^ 8meV/(I//10nm)^, where the 
reduced mass is /i = \{me -1- rrihh) (with the elec¬ 
tron mass rrie = 0.067mo and the heavy-hole mass 
'^hh = O.llmo, and mo is the bare electron mass). 
For vdW heterostructures the effective exciton mass 
is approximately mo [27], so the kinetic energy scale 
is estimated by fi^/^moT^ 150meV/(I//lnm)^. 

We find the bound state of the two IXs by nu¬ 
merically solving Eq. (2). An exemplary solution of 
the probability density |'0(p)p in the lowest energy 
bound state is shown in Fig. 2(b). Below we dis¬ 
cuss the dependence of the binding energy of IXm 
on a range of (realistic) geometrical parameters for 
the two studied exemplary material systems (GaAs 
and vdW). For the symmetric configuration di = ^ 2 , 
the calculated binding energy as a function of the 
dipoles separation L is shown in Fig. 3(a,b) for three 
different dipole sizes. The binding energy increases 
(in absolute value) with decreasing inter-layers sep¬ 
aration L and larger dipole size d, as the attraction 
between two aligned dipoles scales like dP/L^. In ad¬ 
dition, the insets of Fig. 3(a,b) show the full width at 
half maximum (FWHM) of the probability density, 
which is a measure for the molecules’ Bohr radius. 
One can see that the FWHM is of the same order 
of L and increases with it. The calculation is per¬ 
formed in a limited range of inter-bilayer separation 
L, since for L values that are too small the binding 
energy tends to exceed that of the unbound IX it¬ 
self, thus making our treatment inaccurate. Further¬ 
more, from the practical side, making L too small 
might increase the electron inter-bilayer tunneling 
rate. 

Next we study the effect of asymmetry between 
the two bilayers. Asymmetric layers are of practical 
interest and importance since different layer thick¬ 
nesses leads to different exciton energies so that the 
excitons in the two layers can be optically addressed 
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FIG. 3. (a,b) Binding energy of IXm molecule as a func¬ 
tion of the dipoles separation L in GaAs and vdW mate¬ 
rials. (a) di = (i 2 = 120A (blue, solid), di = 0^2 = 140A 
(red, dash-dot) and di = 0^2 = I 6 OA (yellow, dashed), 
(b) di — d 2 — 1.2nm (blue, solid), di — d 2 — 1.4nm 
(red, dash-dot) and di = 0^2 = 1 . 6 nm (yellow, dashed). 
Insets: full width at half maximum (FWHM) of the prob¬ 
ability density |' 0 (p)|^, measured in units of L. (c,d) 
Binding energy of IXm molecule formed in GaAs and 
vdW materials as a function of the structure asymme¬ 
try, (di — d 2 )/ 2 L at fixed L = 300A (c) and L = 3nm 
(d). In both panels (di -hd 2 )/ 2 L was fixed to 6/15 (blue, 
solid), 7/15 (red, dash-dot) and 8/15 (yellow, dashed). 


and resolved independently. Fig. 3(c,d) shows the 
calculated binding energy as a function of the differ¬ 
ence di — (^ 2 , for a few fixed values of di -1- d 2 and L. 
The dependence of binding energy on di ±d 2 at fixed 
L is qualitatively understood as follows: the overall 
increase in binding energy with di + d 2 is explained 
by the growing attraction between the electron and 
the hole in the upper and the lower bilayer respec¬ 
tively, given explicitly by the last term of Eq. (1). 
The decrease in binding energy with di — d 2 at fixed 
di +d 2 and L is because the attraction part above re¬ 
mains unaffected by di — d 2 , while the repulsion part 
contained in the first two terms of Eq. (1) grows. It 
follows that symmetric structures yields the maxi- 


3 




















FIG. 4. (a) A system of three IX’s, one resides in the up¬ 
per DQW and the two others are in the lower DQW, with 
lateral distances of ri 2 , ris, r 23 respectively. The optimal 
configuration minimizing the repulsion for given values 
of ri 2 and ris requires that the three excitons reside on 
the same straight line (when viewed from the top), such 
that r 23 = ri 2 + ri 3 . (b) The interaction potential of 
three dipolar excitons Usir 12 , r 13 , r 23 = '^12 + r 13 ). The 
minimum is reached for ri 2 (i 3 ) ^ 00 and ri 3 ( 23 ) = 0. 

mum binding and are in principle favorable, how¬ 
ever practical considerations might require a slightly 
asymmetric structure. 

Next, we consider the possibility of formation of 
bound complexes that contain more than two IXs. 
Consider the situation depicted in Fig. 4(a), where 
one IX resides in the upper bilayer at a lateral loca¬ 
tion ri , and two other IXs reside in the lower bilayer 
at lateral locations r 2 and respectively. Denoting 
Gj = we can describe the total electrostatic 

interaction energy by 

U3{ri2, ri3, r23) = Ud{ri2) + Ud{ri3) 

^ /q_ 1 \ 

e \r23 y/rls+dlj ’ 

The first two terms are the interaction between the 
two IXs at different bilayers, Eq. ( 1 ), while the third 
term is the same-bilayer dipole repulsion. This last 
term is minimized for fixed ri 2 and ris when all 


three excitons are arranged along a single straight 
line when the system is viewed from above, i.e. 
when ri — r 2 || ri — rs. In this case we have 
^23 = '^12 + ^ 13 - In particular, for a symmetric con¬ 
figuration ri 2 = ri 3 = ^ 23/2 = r, the 3-body po¬ 
tential t/ 3 (r, r, 2 r) reaches a minimum U^rnin < 0 at 
some value of r = r^. This might indicate that a 
bound complex can exist. However, as can be seen in 
Fig. 4, Usrnin < 0 is uot a global potential minimum 
but rather a saddle point, and the true minimum 
is achieved when one of the IXs escapes to infinity. 
Therefore, a complex of three IXs is impossible. On 
the contrary, stacking more than two bilayer struc¬ 
tures vertically, may result in more complex bound 
states of more than two IXs. Since, the same ar¬ 
gument applies to complexes with more than three 
IXs, we conclude that in double bilayer structures 
bound states of more than two excitons do not ex¬ 
ist. Therefore only unbound IXs and IXm should be 
considered here. 

Finally, we estimate the relative fraction of IXm 
under realistic quasi-equilibrium conditions, in the 
dilute exciton gas limit where interactions are neg¬ 
ligible and a Boltzmann distribution applies. Under 
these assumptions it is possible to derive a Saha-like 
equation (see appendix): 

^ 27r?i^ ^E/ksT 
nin2 fiksT 

where ni,n 2 are the unbound excitons densities in 
the two separate bilayers respectively, um is the 
molecules density and E is their binding energy (in 
absolute value). In Fig. 5 we plot the fraction of 
molecules out of the total exciton density as a func¬ 
tion of temperature for the case ni = n 2 . There is 
a significant IXm population at temperatures and 
densities readily achievable experimentally in both 
material systems. Therefore, a clear experimental 
signature of molecules formation is expected in the 
form of a new emission line in the photolumines¬ 
cence spectrum: the recombination of one exciton 
in a molecule will result in a red shifted photolu¬ 
minescence energy relative to unbound IXs due to 
the finite binding energy of the molecule. Owing to 
the residual kinetic energy of the remaining exciton, 
this red shift might be larger than just the molecular 
binding energy. As the calculated values of E can 
be larger than the IX linewidth, this emission should 
be easily detectable. 

In conclusion. We found that vertically coupled 
dipolar molecules can form in vertically coupled bi¬ 
layer QWs, with a typical binding energy of few 
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FIG. 5. The relative fraction of IXm density as a func¬ 
tion of temperature for different total IX densities rit = 
Til + tim = ^2 + riM, in the specific case where ni — n 2 . 
Solid lines: a GaAs-based system with E — O.lmel/. 
Dashed lines: a vdW structure with E — 2meV. 


Kelvins for GaAs-based systems and up to tens 
of Kelvins in vdW based systems, and that such 
molecular species should have a significant popula¬ 
tion within a realistic range of experimental param¬ 
eters, allowing a clear spectroscopic signature. Such 
systems can open up opportunities for establishing 
quantum fluids with more complex interactions, as 
well as means for realizing dipolar drag and dipolar 
chains between two or more stacked layers of dipoles. 
As the observation of single IX in a trap was recently 
reported [28] in a single GaAs bilayer structure, ob¬ 
serving single trapped IXm molecules in a similar 
trap geometries but with a vertically coupled bilayer 
structure could be the next step to realize complex 
exciton manipulations and study dipolar molecular 
dynamics in a well controlled system. 
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Appendix: equilibrium density of IXm 


In the case of low density both exciton and exciton molecule gases are non-degenerate, and multi-particle 
interactions are small [9]. The free energy of exciton gases in the first and second bilayers are 


Fi(ni,T) = — AniTln 
F2{n2,T) = —An2T\n 


2emT 

ni 

2emT 
ttTi^ n2 


(Ala) 

(Alb) 


where A is the cloud area and ni,n 2 are the exciton densities in the different layers respectively (here e is 
the natural logarithm base and a unit system with kB = 1 is utilized). Compared to a single-component 
gas, a factor of 4 appears in the argument of the logarithm due to 4 possible spin states of excitons. 
Similarly, it is possible to write the free energy of exciton molecular gas with density nyf- 


Tm(^M, T) =In--2 - , (A2) 

Trn riM 

where E is the molecular binding energy. Note that there are 16 possible spin state of exciton molecules and 
that the mass of a molecule is 2m. 

Excitons in both layers can be pumped separately and their densities are independent. The total density 
of excitons in each layer is then 


nit = ni + riM , n2t = n2 E um • (A3) 

Given nu and n 2 t, the condition of equilibrium is a minimum of the total free energy, 

F = Fi + F2 + Fm , (A4) 

that is the chemical potential of the molecular gas equals the sum of chemical potentials of the exciton gases: 


where 


— Mm 5 


1 dFi 2mT 

Ml = M — 9 — : 

A dni TrM^ni 

1 dF 2 2mT 

A dn2 ttTl 712 

1 OFm ^ rr, 16mT 

l^M = ^ A = -F - T In —32 - 

ttM um 


A duM 


Substitution of Eq.(A6) in Eq.(A5) gives 


nxM = 


Aira 

mT 


e^/'^nxinx2 ■ 


(A5) 

(A6a) 

(A6b) 

(A6c) 

(A7) 


Substituting the reduced mass // = m/2 yields Eq.(4) of the main text. Elimination of ni and n 2 with the 
help of Eq.(A3) leads to an equation for tim' 


2 / N niT 

- (^it + n2t)nM - 

AttTi 


e + nitn 2 t = 0 . 


(A8) 
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For example, exciton concentrations of nu = n 2 t = lO^cm ^ meet all of the above conditions and give 
exciton molecule concentration of tim ~ 0.5 x 10^cm“^. 
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